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Abstract 

The Lorentz-invariant S-matrix elements in interacting quantum field theory 
(QFT) are used to represent the QFT state by a Lorentz-invariant many-time 
wave function. Such a wave function can be used to describe inelastic scattering 
processes (involving particle creation and destruction) by Bohmian particle 
trajectories satisfying relativistic-covariant equations of motion. 

"^ ! PACS Numbers: ll.10.-z, 11.80.-m, 03.65.Ta 
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2 : 1 Introduction 



The Bohmian formulation of quantum mechanics (QM) [U EJ US IU G3 E] allows to 
visualize quantum processes in terms of continuous deterministic pointlike-particle 
trajectories guided by wave functions. In this formulation, all quantum uncertainties 
emerge from ignorance of the actual initial particle positions. Such a formulation not 
only offers a possible answer to deep conceptual and interpretational puzzles of QM, 
but in some cases provides also a powerful practically useful computational tool [3j [7] . 
A challenge for the Bohmian formulation of quantum theory is to reconcile con- 
tinuous particle trajectories with phenomena of particle creation and destruction in 
quantum field theory (QFT). One possibility is to explicitly break the rule of continu- 
ous deterministic evolution, by adding an additional equation that specifies stochastic 
breaking of the trajectories |S1 19]. Another possibility is to introduce an additional 
continuously and deterministically evolving hidden variable that specifies effect ivity 
of each particle trajectory (TUJ [IT]. However, both possibilities seem rather artificial 
and contrived. In addition, the explicit constructions in [HI EJ EH E] do not obey 
Lorentz invar iance. 



Recently, a much simpler approach has been introduced [12] , in which the particle 
trajectories are continuous, but the appearance of particle creation and destruction 
results from quantum theory of measurements describing entanglement with particle 
detectors. The QFT states are represented by many-time wave functions (which, 
as a byproduct, offers a new resolution of the black- hole information paradox [13]), 
implying that the corresponding Bohmian equations of motion for particle trajectories 
are relativistic-covariant. The equations determining the many-time wave functions 
are explicitly Lorentz-covariant for free fields. In the interacting case, however, the 
many-time evolution of wave functions in [12] is described by the time-evolution 
operators U(t) in the interaction picture, which lack manifest Lorentz invariance for 
interacting QFT. 

In this paper we further develop some of the ideas introduced in [12] . In particular, 
we find a new method for calculation of the many-time wave function in interacting 
QFT, which turns out to be (i) simpler than that in [T2] and (ii) manifestly Lorentz- 
invariant. In addition, we propose a new guiding equation for Bohmian trajectories 
which turns out to be much simpler than that in [12] . We also clarify some conceptual 
issues in a somewhat different (and hopefully more illuminating) way than in [T2] , 
and discuss some limitations of the present approach that represent a challenge for 
the future research. 

The paper is organized as follows. In Sec. [2] we present a brief review of the 
Bohmian formulation of relativistic QM, as well as a general qualitative description of 
particle creation and destruction, with emphasis on conceptual (rather than technical) 
aspects. In Sec. [3] we describe in detail how states in free and interacting QFT can 
be represented by Lorentz-invariant many-time wave functions. In Sec. H] we give a 
probabilistic interpretation of these wave functions and develop the corresponding 
Bohmian interpretation in terms of particle trajectories. The conclusions are drawn 
in Sec. |5j 

In the paper we use units h — c — 1 and the metric signature (+, — , — , — ). 

2 Conceptual preliminaries 

2.1 Outline of relativistic Bohmian mechanics 

Since Bohmian mechanics (BM) is a nonlocal theory, it is frequently objected that it 
is not compatible with relativity. Yet, as reviewed in [T3], all objections of that type 
can be circumvented. Relativistic-covariant nonlocal Bohmian equations of motion 
for particle trajectories have been introduced in [15] and further studied in [HI [TO 
IT8] . Using the idea that, in relativistic QM, space probability density should be 
generalized to spacetime probability density (see, e.g., [191 ED]) , a relativistic-invariant 
probabilistic interpretation (associated with relativistic particle trajectories) has been 
introduced in [201 IZI] ■ As indicated in [12] and generally shown in [T3] , this makes the 
measurable statistical predictions of relativistic BM compatible with all measurable 
statistical predictions of the "ordinary" purely probabilistic interpretation of QM. 
Let us present a brief overview of the main ideas of relativistic BM. Denoting a 



point in spacetime by x = {x M }, an n-particle quantum state is represented by a many- 
time wave function ip(xi, . . . ,x n ). This wave function lives in the 4n- dimensional 
configuration space with coordinates x%, a — 1, ... ,n, referred to as relativistic con- 
figuration space. This is a many-time wave function [22J, because each particle has 
its own time-coordinate x Q a . The quantity i[)*ip is interpreted as probability density 
in the relativistic configuration space, in the sense that the infinitesimal probability 
dP of finding n particles in an infinitesimal 4n- dimensional volume around the points 
Xx,...,x n is given by 

dP = ip*(xi, ..., x n )ip(xi, . . . , x n )d 4 xi • ■ ■ d A x n . (1) 

The wave function satisfies the n-particle Klein-Gordon equation 

£[(9^ + 77^ = 0, (2) 

o=l 

which implies that the n-particle Klein- Gordon current 

J"a = \rk^ (3) 

is conserved: 

n 

E^^ M = o. (4) 

a=l 

The Bohmian particle trajectories are integral curves of the vector field V£(xi, . . . ,x r 
calculated from the wave function as 

v ' = & < 5 > 

In the parameterized form, the integral curves are represented by functions X£(s) 
satisfying 

—^ r 1 = V»(X 1 (s),...,X n (s)). (6) 

From (jlj) and the fact that ip and does not have an explicit dependence on s, one 
finds the equivariance equation 

^ + EW«') = ». (7) 

os o=l 

Formally, Eq. (j7]) shows that ([6]) is compatible with ([1]) in the following sense: If a 
statistical ensemble of particles has the probabilistic distribution ([1]) for some initial 
s, then it has the probabilistic distribution (JT]) for any s. 

To understand the physical meaning of the formal equivariance equation ([7j), one 
needs to understand the physical meaning of the parameter s. For that purpose, it 
is useful to exploit the analogy with nonrelativistic Newtonian mechanics [23] ■ For 
nonrelativistic particle systems with conserved energy, the forces do not have an 
explicit dependence on time t. The only quantities that have a dependence on t are 



particle trajectories X^(t), i = 1,2,3. Thus, the parameter t has a physical meaning 
only along trajectories; time without trajectories does not exist! In this sense, t is only 
an auxiliary parameter that serves to parameterize the trajectories in 3-dimensional 
space, not a fundamental physical quantity by its own. Yet, a "clock" can measure 
time indirectly. Namely, a "clock" is nothing but a physical process described by a 
function X l a {t) periodic in t. One actually observes the number of periods, and then 
interprets it as a measure of elapsed time. 

The theory of relativity revolutionized the concept of time by replacing the pa- 
rameter t with a coordinate x° treated as a 4th dimension not much different from 
other 3 space dimensions. Yet, it does not mean that an auxiliary Newton-like time 
parameter is completely eliminated from relativistic mechanics. Such a parameter 
can still be introduced to parameterize relativistic spacetime particle trajectories in 
a manifestly covariant manner. This parameter, denoted by s, can be identified with 
a generalized proper time defined along particle trajectories of many-particle systems 
|23j . The parameter s can even be measured indirectly by a "clock" corresponding 
to a physical process periodic in s, in complete analogy with measurement of t in 
nonrelativistic mechanics. As discussed in more detail in [23], this makes the pa- 
rameter s appearing in fl6]) a physical quantity, very much analogous to Newton time 
t. With this physical insight, the relativistic-covariant equation ([7]) is to be inter- 
preted as physical probability conservation during the evolution parameterized by 
the evolution-parameter s. 

Now let us sketch how all statistical predictions of the purely probabilistic inter- 
pretation of QM can be reproduced from relativistic BM (for more details see [H]). 
Let a physical system be described by a wave function 

*l>( x ) = J2°bM x ), ( 8 ) 

b 

where ipb{x) are eigenstates of some hermitian operator B on the Hilbert space of func- 
tions of x, normalized in a large but finite 4-dimensional box such that / d 4 xipl(x)ifji,(x) 
1. The purely probabilistic interpretation asserts that |q,| 2 is the probability that the 
observable B will take the value b. To see how BM reproduces this assertion, one 
needs to take into account the entanglement with the measuring apparatus. This 
leads to a wave function [TJ] 

^fo y) = J2 c b A( x ) E b(y), (9) 

b 

where Eb(y) are detector wave functions normalized in the relativistic configuration 
space, with a negligible overlap in that space. Since §6§ is compatible with ([I]), 
the probability that detector particles will have a relativistic configuration Y in the 
support of Eb(y) is equal to |q,| 2 [H]. This shows that the measurable statistical pre- 
dictions of relativistic BM coincide with those of purely probabilistic interpretation. 
In particular, if ipt,{ x ) are (approximate) eigenstates of the space-position opera- 
tors x 1 , x 2 , x 3 at time x°, then the measurement procedure above (approximately) 
reproduces the usual space probability distribution dP{z) = ip*(x)il)(x)d 3 x. 



2.2 The general mechanism of particle creation and destruc- 
tion 

In Sec. 12.11 we have studied relativistic QM with a fixed number n of particles. Now 
we give a qualitative description of the general mechanism of particle creation and 
destruction in QFT. Schematically, it can be described in 3 steps. 

In the first step, the initial state Initial) with a definite initial number of particles 
^initial suffers a unitary deterministic evolution in interacting QFT 

^initial) ->^2c n \n), (10) 

n 

where the final state is a superposition of states \n) with different numbers of particles. 
In the second step, the quantum state above interacts with the environment (e.g., 
a particle detector), which causes a unitary evolution that creates entanglement with 
environment states \E n ): 



J2 c n\n) 



^initial) -> ^C n \n)\E n ). (11) 



Here state \E n ) can be thought of as a macroscopic state describing a detector in 
a state of saying that n particles are detected. Since different macroscopic states 
are macroscopically distinguishable, the corresponding wave functions E n (y) have a 
negligible overlap: 

E n (y)E n ,(y)~0 for n^n'. (12) 

In the third step, one needs a mechanism that will pick up only one term in 
the superposition (TlTj) . Conventionally, it is usually described by the wave-function 
"collapse" . The role of the Bohmian formulation is to replace this ad hoc collapse 
with a mathematically better defined physical process. The wave function depending 
on y guides the detector particles with a trajectory Y(s). Due to (TT21 . the particles 
enter only one channel £ , riflnal (?/) among many channels E n (y) in (TTTj) . This makes all 
other channels empty, which for all practical purposes is effectively the same as if the 
state exhibited a collapse 



E 



c n \n)\E n ) -> |nfi nal )|£ nflnal ). (13) 



(For more details see also [12].) 

The effect of these 3 steps can be summarized as a transition 

Initial) | -^initial) ~ > | ""-final ) | -En flnal ) , (14) 

which typically involves the destruction of some initial particles and the creation of 
some new (final) ones. 



3 QFT states represented by Lorentz- invariant wave 
functions 

3.1 Momentum-eigenstates with fixed number of particles 

Free QFT is usually formulated in terms of n-particle states \ki,...,k n ) with on- 
shell 4-momenta k a , a = l,...,n. (For simplicity, we consider particles without 
spin, but spin can also be included as in [12J.) Introducing the condensed notation 
/?") = (fei, . . . , k n ), we denote these states as 

|n,^) = |A;i,...,A; n }. (15) 

The state \n, k^) can also be represented by an n-point wave function [T2] 

y n i^(n){Xx, . . . , x n ) = \Xi, ■ ■ . , x n \ fci, . . . , K n ) 

= b n S {xi _ Xn} e- ik ^---e- lk - x ", (16) 

where b n is a normalization factor and Si Xlr .. tX \ denotes symmetrization over all 
x\, . . . , x n . The normalization factor b n is chosen such that [12] 



Vx^ nrkin) (x^)\ 2 = l, (17) 

where x^ = {x\, . . . , x n ) and 

Vx (n) = d 4 x 1 ---d A x n . (18) 

3.2 General states in free QFT 

A general free QFT state has an expansion of the form 

oo 

l^) = co|o) + EE^(,)K^ (n) ), (i9) 

n=l £(n) 

where |0) is the vacuum. We assume that the state is normalized to unity, in the 
sense that 

oo 

W> = |c | 2 + EEk,^)l 2 = i- (20) 

n=l %( n ) 

Since (fT§|) involves a superposition of states with various numbers n of particles, 
it is convenient to slightly modify the notation of Sec. 13.11 To distinguish particle 
positions x a coming from sectors of different n, instead of x^ = (xx, . . . , x n ) we write 

X = \X n ^i, • • • , Xn^nJ. v^-*-/ 

Then the collection of particle positions from sectors of all possible n is denoted as 

x = (f (1) ,f (2) ,...) = (x lt i,x 2 ,i,x 2 ,2,---)- (22) 
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The state (TT9|) can be represented by a many-component wave function 



ty(x) 



^i(f«) 
^ 2 (£ (2) ) 



\ 



(23) 



/ 



Here 



*r 



v c °' 



^jf (n) : 






where 



/oo 
ra=l 



are volumes of the corresponding relativistic configuration spaces, and 



1pn(x n ) — Z^ C n ) fe(n)V'n,^(™) ' 



X 



(n^ 



(24) 
(25) 

(26) 
(27) 



are n-particle wave packets with ^ n t(n){x^ n ') given by (TT6l) . The tilde above wave 
functions denotes wave functions which are not necessarily normalized to unity. The 
normalization factors in (T24|) - (T25|) are chosen such that [T2] 



£>£|^nl 2 



l c o| 2 , 






|2 



where 



Vx = \[Vx {n) . 



(2t 



(29) 



71 = 1 



This provides that the total wave function (1231) is normalized to unity, in the sense 
that 

fvx¥{x)^{x) = fvx\^ \ 2 + ]T fvx\^ n \ 2 = 1. (30) 

n=l 

Thus we see that the many-time wave function in free QFT is uniquely determined 
by a set of expansion coefficients Co, c n ^ (n) . 

3.3 Scattering wave function 

In interacting QFT, particles may be created and destructed. The information on 
dynamics of this creation and destruction in encoded in the scattering matrix (shortly, 
S-matrix) with matrix elements in the momentum space (see, e.g., [24" t [25 | [26]) 



(h, ..., k n \S\ Pl , ..., Pm ) = S(n, k^\m,^). 



(31) 



These matrix elements are Lorentz- invariant. In general, the initial state before scat- 
tering is a superposition of the form (TT9l) with expansion coefficients cfj 1 , c m ~,y Then, 
the final state after scattering is also a superposition of the form f fl9|) . but with dif- 
ferent expansion coefficients Cp ut , c ou i {n) given by 

c 3(«) = E E S(n, fcW |m, j*"0) cj^ m) . (32) 

(Eq. ( 132|) can be thought of as including the vacuum terms with n = and m = 0, 
but they are not written explicitly because the S-matrix elements involving vacuum 
are usually trivial.) Thus, the total wave function can be written as 



%f(x) 



-^ in (£) for xeTZ m , 

F- (33) 

V ^ out (f) foTxeTZ out , 



V 



where 7£ m and TZ ont are in and out regions of the relativistic configuration space 
with volumes V in and y out , respectively. The wave functions ty m (x) and \l/ out (x) are 
defined as in Sec. 13.21 with coefficients cJJ 1 , c m - (n) and Cq u1 , c on % n) , respectively. The 

normalization factors in ( 1331) are chosen such that fDxW(x)^(x) = 1. 

Let us make a few comments on validity of the approximation ( 1331) . Strictly 
speaking, the S-matrix elements (I3TJ) refer only to in states \p\, . . . ,p m ) at x° a — > — oo 
and out states \ki,...,k m ) at x° a — > oo. Nevertheless, they also represent a good 
approximation for large but finite in and out times. After all, the predictions obtained 
from these matrix elements are in excellent agreement with experiments, and no 
experiments are performed at infinity. In fact, the S-matrix elements represent a 
good approximation wherever particles can be well approximated by free particles. 
But this means that the approximation ( 133]) is good almost everywhere, except in a 
very small region of spacetime at which the collision of localized wave packets actually 
happens. For relativistic collisions, the size of this small region of collision is typically 
of the order 1/E, where E is a typical energy of the colliding particles. (For decay 
processes, E is the decay width, which is large for short-living particles.) Thus, we can 
conclude that, in a typical situation, ( 133|) is a good approximation almost everywhere, 
i.e., that lZ m U lZ out covers almost the whole relativistic configuration space. 

Let us compare it with the results of [12]. In principle, the method developed in 
[T2] allows to calculate \l/(x) everywhere, but is not manifestly Lorentz- invariant. By 
contrast, the method developed here is manifestly Lorentz- invariant, but is not valid 
everywhere. It seems that it should be possible to develop a method that is both 
valid everywhere and manifestly Lorentz- invariant, but we leave the development of 
such a method as a program for the future research. 

4 Probability and particle trajectories in QFT 

In this section we develop the physical interpretation of the wave function \I/(x), as a 
natural generalization of the results in Sec. 12.11 The probabilistic interpretation (pQ) 



generalizes to 

VP = ^(x)V(x)Vx. (34) 

Due to ( |23|) . the probability density W^ decomposes as 

oo 

¥(x)V(x) = |* ol 2 + E \^n(x^)\ 2 . (35) 

n=l 

Each n-particle wave function ^f n (x^) satisfies the n-particle Klein-Gordon equation 
of the form of ([2]) . 

For further analysis, it is convenient to introduce a condensed label A = (n,a n ), 
such that (1221) can be written as 



x = {x A } = (x 1 ,x 2} x 3 ,...). (36) 

With this notation, we introduce the current 

J1 = ^VL*, (37) 

which generalizes ([3]). Due to fl35|) . it can be decomposed into a collection of n-particle 
currents, each of which is conserved due to the n-particle Klein-Gordon equation. This 
implies that ( 13"T|) is also conserved 

oo 

E d ^ J A = 0, (38) 

A=l 

which implies a generalization of ([7]) 

^|^ + E M^Wjf) = 0, (39) 

as A=l 

where 

« = il < 40 > 

generalizes (jSJ). The Bohmian particle trajectories are given by the generalization of 
© 

^M = ^(Xi(s),* 2 (*),XsOO, • • •), (41) 

which are compatible with ( 134")) due to the equivariance equation (1391) . 

It is straightforward to apply this to the scattering wave function (|33|) . which, to- 
gether with physical insights from Sec. 12.21 provides a Bohmian description of inelastic 
scattering processes involving particle creation and destruction. 



5 Conclusion 

The Bohmian formulation of quantum theory describes all quantum processes in terms 
of continuous deterministic particle trajectories guided by wave functions. The results 
of the present paper show that it can be formulated in a form which (i) obeys manifest 
Lorentz invariance and (ii) includes a description of particle creation and destruction 
in QFT. (For simplicity, in this paper we have discussed only particles without spin, 
but, by using the results presented in p2], the generalization to particles with spin is 
relatively simple.) These results reinforce the view that the Bohmian formulation is 
a viable formulation of quantum theory. 
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